Potential Capabilities of Lunar Laser Ranging for 
Geodesy and Relativity 



Jiirgen Miiller 

Institut fur Erdmessung, University of Hannover, Schneiderberg 50, 30167 Hannover, Germany 

mueller @ ife . uni-hanno ver. de 

James G. Williams, Slava G. Turyshev 

Jet Propulsion Laboratory, 4800 Oak Grove Drive, Pasadena, CA 91 109, USA 
Peter J. Shelus 

University of Texas at Austin, Center for Space Research, 3925 W. Braker Lane, Austin, TX 78759, USA 



in 
o 
o 

(N 
o 

Q 

m 

> 
On 

o 

On 
O 

in 
o 

"o 

7' 
u 



Abstract. Lunar Laser Ranging (LLR), which has 
been carried out for more than 35 years, is used to 
determine many parameters within the Earth-Moon 
system. This includes coordinates of terrestrial rang- 
ing stations and that of lunar retro-reflectors, as well 
as lunar orbit, gravity field, and its tidal accelera- 
tion. LLR data analysis also performs a number of 
gravitational physics experiments such as test of the 
equivalence principle, search for time variation of the 
gravitational constant, and determines value of sev- 
eral metric gravity parameters. These gravitational 
physics parameters cause both secular and periodic 
effects on the lunar orbit that are detectable with 
LLR. Furthermore, LLR contributes to the determi- 
nation of Earth orientation parameters (EOP) such 
as nutation, precession (including relativistic preces- 
sion), polar motion, and UT1. The corresponding 
LLR EOP series is three decades long. LLR can 
be used for the realization of both the terrestrial and 
selenocentric reference frames. The realization of a 
dynamically defined inertial reference frame, in con- 
trast to the kinematically realized frame of VLBI, of- 
fers new possibilities for mutual cross-checking and 
confirmation. Finally, LLR also investigates the pro- 
cesses related to the Moon's interior dynamics. 

Here, we review the LLR technique focusing on 
its impact on Geodesy and Relativity. We discuss the 
modern observational accuracy and the level of ex- 
isting LLR modeling. We present the near-term ob- 
jectives and emphasize improvements needed to fully 
utilize the scientific potential of LLR. 

Keywords. Lunar Laser Ranging, Relativity, Earth- 
Moon dynamics 



1 Motivation 

Being one of the first space geodetic techniques, lu- 
nar laser ranging (LLR) has routinely provided ob- 



servations for more than 35 years. The LLR data are 
collected as normal points, i.e. the combination of 
round trip light times of lunar returns obtained over 
a short time span of 10 to 20 minutes. Out of « 10 19 
photons sent per pulse by the transmit ter, less than 1 
is stat istically detected at the receiver Willi ams et al. I 
this is because of the combination of sev- 
eral factors, namely energy loss (i.e. 1/R 4 law), at- 
mospherical extinction and geometric reasons (rather 
small telescope apertures and reflector areas). More- 
over, the detection of real lunar returns is rather diffi- 
cult as dedicated data filtering (spatially, temporally 
and spectrally) is required. These conditions are the 
main reason, why only a few observatories world- 
wide are capable of laser ranging to the Moon. 

Observations began shortly after the first Apollo 
11 manned mission to the Moon in 1969 which 
deployed a passive retro-reflector on its surface. 
Two American and two French-built reflector ar- 
rays (transported by Soviet spacecraft) followed un- 
til 1973. Since then over 16,000 LLR measure- 
ments have by now been made of the distance be- 
tween Earth observatories and lunar reflectors. Most 
LLR data have been collected by a site operated by 
the Observatoire de la Cote dAzur (OCA), France. 
The transmitter/receiver used by OCA is a 1 .5m alt- 
az Ritchey-Chretien reflecting telescope. The mount 
and control electronics insure blind tracking on a lu- 
nar feature at the 1 arcsec level for 10 minutes. The 
OCA station uses a neodymium-YAG laser, emitting 
a train of pulses, each with a width of several tens 
of picoseconds). LLR station at the McDonald Ob- 
servatory in Texas, USA is another major provider of 
the LLR data. The McDonald Laser Ranging Station 
(MLRS) is built around a computer-controlled 0.76- 
m x-y mounted Cassegrain/Coude reflecting tele- 
scope and a short pulse, frequency doubled, 532- 
nm, neodymium-YAG laser with appropriate com- 
puter, electronic, meteorological, and timing inter- 
faces. Until 1990, Haleakala laser ranging station 



on the island of Maui (Hawaii, USA) contributed to 
LLR activities with its 0.40 m telescope. Single lu- 
nar returns are available from Orroral laser ranging 
station in Australia (closed 1 November 1998) and 
the Wettzell Laser Ranging System (0.75 m) in Ger- 
many. Other modern stations have demonstrated lu- 
nar capability, e.g. the Matera Laser Ranging Sta- 
tion (0.50 m) in Italy and Hartebeesthook Observa- 
tory (0.762 m) in South Africa. A new site with lunar 
capability is currently being built at the Apache Point 
Observatory (New Mexico, USA) around a 3.5 m 
telescope. This station, called APOLLO, is designed 
for mm accuracy ranging (Williams et al. 2004b). 

Today MLRS and OCA are the only currently op- 
erational LLR sites achieving a typical range pre- 
cision of 18-25 mm. Fig. 1 shows the number of 
LLR normal points per year since 1970. As shown 
in Fig. 2, the range data have not been accumulated 
uniformly; substantial variations in data density exist 
as a function of synodic angle D, these phase angles 
are represented by 36 bins of 10 degree width. In 
Fig. 2, data gaps are seen near new Moon (0 and 360 
degrees) and full Moon (180 degrees) phases. The 
properties of this data distribution are a consequence 
of operational restrictions, such as difficulties to op- 
erate near the bright sun in daylight (i.e. new Moon) 
or of high background solar illumination noise (i.e. 
full Moon). 
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Fig. 1. Lunar observations per year, 1970 - 2005. 

While measurement precision for all model pa- 
rameters benefit from the ever-increasing improve- 
ment in precision of individual range measurements 
(which now is at the few cm level, see also Fig. 3), 
some parameters of scientific interest, such as time 
variation of Newton's coupling parameter G/G or 
precession rate of lunar perigee, particularly benefit 
from the long time period (35 years and growing) of 
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Fig. 2. Data distribution as a function of the synodic angle D. 
range measurements. 
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Fig. 3. Weighted residuals (observed-computed Earth-Moon 
distance) annually averaged. 

In the 1970s LLR was an early space technique 
for determining Earth orientation parameters (EOP). 
EOP data from LLR are computed for those nights 
where sufficient data were available (approximately 
2000 nights over 35 years). Although, the other 
space geodetic techniques (i.e., SLR, VLBI) dom- 
inate since the 80s, today LLR still delivers very 
competitive results, and because of large improve- 
ments in ranging precision (30 cm in 1969 to 2 cm 
today), it now serves as one of the strongest tools 
in the solar system for testing general relativity. 
Moreover, parameters such as the station coordinates 
and velocities contributed to the International Ter- 
restrial Reference Frame ITRF2000, EOP quantities 
were used in combined solutions of the International 
Earth Rotation and Reference Systems Service IERS 



(er = 0.5 mas). 

2 LLR Model 

The existing LLR model has been developed to com- 
pute the LLR observables — the round trip travel 
times of laser pulses between stations on the Earth 
and passive reflectors on the Moon (see e.g. Miiller et 
al. 1996, Miiller and Nordtvedt 1998 or Miiller 2000, 
2001, Miiller and Tesmer 2002, Williams et al. 2005b 
and the references therein). The model is fully rel- 
ativistic and is complete up to first post-Newtonian 
(1/c 2 ) level; it uses the Einstein's general theory of 
relativity - the standard theory of gravity. The mod- 
eling of the relativistic parts is much more challeng- 
ing than, e.g., in SLR, because the relativistic correc- 
tions increase the farther the distance becomes. The 
modeling of the 'classical' parts has been set up ac- 
cording to IERS Conventions (IERS 2003), but it is 
restricted to the 1 cm level. Based upon this model, 
two groups of parameters (170 in total) are deter- 
mined by a weighted least- squares fit of the observa- 
tions. The first group comprised from the so-called 
1 Newtonian' parameters such as 

- geocentric coordinates of three Earth-based 
LLR stations and their velocities; 

- a set of EOPs (luni-solar precession constant, 
nutation coefficients of the 18.6 years period, 
Earth's rotation UT0 and variation of latitude by 
polar motion); 

- selenocentric coordinates of four retro- 
reflectors; 

- rotation of the Moon at one initial epoch (phys- 
ical librations); 

- orbit (position and velocity) of the Moon at this 
epoch; 

- orbit of the Earth-Moon system about the Sun at 
one epoch; 

- mass of the Earth-Moon system times the grav- 
itational constant; 

- the lowest mass multipole moments of the 
Moon; 

- lunar Love number and a rotational energy dis- 
sipation parameter; 

- lag angle indicating the lunar tidal acceleration 
responsible for the increase of the Earth-Moon 
distance (about 3.8 cm/yr), the increase in the 



lunar orbit period and the slowdown of Earth's 
angular velocity. 

The second group of parameters used to perform 
LLR tests of plausible modifications of general the- 
ory of relativity (these parameter values for general 
relativity are given in parentheses): 

- geodetic de Sitter precession f^s of the lunar 
orbit (~ 1.92" /cy); 

- space-curvature parameter 7 (= 1) and non- 
linearity parameter (3 (= 1); 

- time variation of the gravitational coupling pa- 
rameter G/G (= yr _1 ) which is important for 
the unification of the fundamental interactions; 

- strong equivalence principle (EP) parameter, 
which for metric theories is r\ = 4/3 — 3 — 7 
(=0); 

- EP-violating coupling of normal matter to 'dark 
matter' at the galactic center; 

- coupling constant a (= 0) of Yukawa potential 
for the Earth-Moon distance which corresponds 
to a test of Newton's inverse square law; 

- combination of parameters Ci — Co — 1 (= 0) de- 
rived in the Mansouri and Sexl (1977) formal- 
ism indicating a violation of special relativity 
(there: Lorentz contraction parameter (1 = 1/2, 
time dilation parameter Co = —1/2); 

- a\ (= 0) and u<i (= 0) which parameterize 'pre- 
ferred frame' effects in metric gravity. 

Most relativistic effects produce periodic pertur- 
bations of the Earth-Moon range 

n 

Ar BM = ^A 4 cos(^At + $ t ). (1) 

i=l 

Ai, uj i} and fa are the amplitudes, frequencies, and 
phases, respectively, of the various perturbations. 
Some example periods of perturbations important for 
the measurement of various parameters are given in 
Table 1 .Note: the designations should not be used as 
formulae for the computation of the corresponding 
periods, e.g. the period 'sidereal-2-annual' has to be 
calculated as l/(l/27.32 d - 2/365.25 d ) w 32.13 d . 
'secular + emerging periodic' means the changing or- 
bital frequencies induced by G/G are starting to be- 
come better signals than the secular rate of change of 
the Earth-Moon range in LLR. 
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Fig. 4. Sensitivity of LLR with respect to G/G assuming 
AG/G = 8 ■ lO" 13 yr" 1 . 
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Fig. 6. Power spectrum of a possible equivalence principle 
violation assuming A(tog /mj) ~ 10 -13 . 
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Fig. 5. Power spectrum of the effect of G/G in the Earth- 
Moon distance assuming AG/G = 8 ■ 10 -13 yr -1 . 



ativistic phenomena show up with typical periods). 
Other effects, like the asteroids cause similar orbit 
perturbations (i.e. with the same frequencies) but 
with different amplitudes, so that they can be distin- 
guished and separated from the effects investigated. 

Table 1. Typical periods of some relativistic quantities, taken 
from Miiller et al. (1999). 



Parameter 


Typical Periods 


V 


synodic (29 d 12'M4 m 2.9 s ) 


Ci-Co-i 


annual (365.25 d ) 


^fi'galactic 


sidereal (27 d 7 h 43 m 11.5 s ) 




sidereal, annual, sidereal-2-annual, 
anomal. (27 d 13' l 18 m 33.2 s )±annual, synodic 


Q2 


2-sidereal, 2-sidereal-anomal., nodal (6798 d ) 


G/G 


secular + emerging periodic 



Fig. 4 represents the sensitivity of the Earth-Moon 
distance with respect to a possible temporal variation 
of the gravitational constant in the order of 8 • 1CP 13 , 
the present accuracy of that parameter. It seems as 
if perturbations of up to 9 meters are still caused, 
but this range (compared to the ranging accuracy at 
the cm level) can not fully be exploited, because the 
lunar tidal acceleration perturbation is similar. The 
largest periods for G/G are shown in Fig. 5 and for 
the EP-parameter in Fig. 6. Obviously many peri- 
ods are affected simultaneously, because the pertur- 
bations, even if caused by a single beat period only 
(e.g. the synodic month for r/), change the whole lu- 
nar orbit (and rotation) and therefore excite further 
frequencies. Nevertheless these properties can be 
used to identify and separate the different effects and 
to determine corresponding parameters (note that rel- 



3 Results 

The global adjustment of the model by least-squares- 
fit procedures gives improved values for the esti- 
mated parameters and their formal standard errors, 
while consideration of parameter correlations ob- 
tained from the covariance analysis and of model 
limitations lead to more 'realistic' errors. Incom- 
pletely modeled solid Earth tides, ocean loading or 
geocenter motion, and uncertainties in values of fixed 
model parameters have to be considered in those es- 
timations. For the temporal variation of the gravita- 
tional constant, G/G = (6 ± 8) • 10 -13 has been ob- 
tained, where the formal standard deviation has been 
scaled by a factor 3 to yield the given value. This 
parameter benefits most from the long time span of 
LLR data and has experienced the biggest improve- 



ment over the past years (cf. Miiller et al. 1999). 
In contrast, the EP-parameter 77 (= (6 ± 7) • 1CT 4 ) 
benefits most from highest accuracy over a sufficient 
long time span (e.g. one year) and a good data cover- 
age over the synodic month, as far as possible. Its im- 
provement was not so big, as the LLR RMS residuals 
increased a little bit in the past years, compare Fig. 3. 
The reason for that increase is not completely under- 
stood and has to be investigated further. In combi- 
nation with the recent value of the space-curvature 
parameter 7 C assini (7 - 1 = (2.1 ± 2.3) • 1(T 5 ) de- 
rived from Doppler measurements to the Cassini 
spacecraft (Bertotti et al. 2003), the non-linearity pa- 
rameter (3 can be determined by applying the re- 
lationship 77 = 4/3 — 3 — 7cassini- One obtains 
/3-l = (1.5 ± 1.8) -lO" 4 (note that using the EP test 
to determine parameters 77 and (3 assumes that there 
is no composition-induced EP violation). 

Final results for all relativistic parameters ob- 
tained from the IfE (Institute fur Erdmessung) anal- 
ysis are shown in Table 2. The realistic errors are 
comparable with those obtained in other recent in- 
vestigations, e.g. at JPL (see Williams et al. 1996, 
2004a, 2004b, 2005b). 

Table 2. Determined values for the relativistic quantities and 
their realistic errors. 



Parameter 


Results 


diff. geod. prec. ^gp - ^dcSit ["/cy] 


(6 ± 10) ■ 10" a 


metric par. 7 — 1 


(4 ± 5) ■ icr :i 


metric par. /3 — 1 (direct) 
and from rj = 4/3 - 3 - 7Ca SS ini 


(-2 ±4) • 
(1.5 ± 1.8) ■ 10" 4 


equiv. principle par. -q 


(6 ± 7) • icr 4 


time var. grav. const. G/G [yr _1 ] 


(6 ±8) • 10- 13 


Yuk. coupl. const. a A=4 . 10 5 km 


(3 ±2) • 10- 11 


spec, relativi. C,i — Co — 1 


(-5 ± 12) ■ 10" & 


infl. of dark matter <5g ga i act i c [erf] 


(4 ±4) • 10- 14 


'preferred frame' effect «i 


(-7 ± 9) • 10- b 


'preferred frame' effect 02 


(1.8 ± 2.5) ■ 10- 5 



4 Further Applications 

In addition to the relativistic phenomena mentioned 
above, more effects related to lunar physics, geo- 
sciences, and geodesy can be investigated. The fol- 
lowing items are of special interest: 

1. Celestial reference frame: A dynamical real- 
ization of the International Celestial Reference 
System (ICRS) by the lunar orbit is obtained 
(cr = 0.001") from LLR data. This can be com- 
pared and analyzed with respect to the kinemati- 



cal ICRS from VLBI. Here, the very good long- 
term stability of the orbit is of great advantage. 

2. Terrestrial reference frame: The results for the 
station coordinates and velocities, which are es- 
timated simultaneously in the standard solution, 
contribute to the realization of the international 
terrestrial reference frame, e.g. to the last one, 
the ITRF2000. 

3. Earth rotation: LLR contributes, among oth- 
ers, to the determination of long-term nutation 
parameters, where again the stable, highly ac- 
curate orbit and the lack of non-conservative 
forces from atmosphere (which affect satellite 
orbits substantially) is very convenient. Addi- 
tionally UT0 and VOL (variation of latitude) 
values are computed (e.g. Dickey et al. 1985), 
which stabilize the combined EOP series, es- 
pecially in the 1970s when no good data from 
other space geodetic techniques were available. 
The precession rate is another example in this 
respect. The present accuracy of the long- 
term nutation coefficients and precession rate 
fits well with the VLBI solutions (within the 
present error bars), see (e.g. Williams et al. 
2005a,c) 

4. Relativity: In addition to the use of LLR in the 
more 'classical' geodetic areas, the dedicated 
investigation of Einstein's theory of relativity 
is of special interest. With an improved accu- 
racy the investigation of further effects (e.g. the 
Lense-Thirring precession) or those of alterna- 
tive theories become possible. 

5. Lunar physics: By the determination of the li- 
bration angles of the Moon, LLR gives access 
to underlying processes affecting lunar rotation 
(e.g. Moon's core, dissipation), cf. Williams et 
al. (2005a). A better distribution of the retro- 
reflectors on the Moon (see Fig. 7) would be 
very helpful. 

6. Selenocentric reference frame: The determina- 
tion of a selenocentric reference frame, the com- 
bination with high-resolution images and the es- 
tablishment of a better geodetic network on the 
Moon is a further big item, which then allows 
accurate lunar mapping. The LLR frame is used 
as reference in many application, e.g. to derive 
gravity models of the Moon (e.g. Konopliv et 
al. 2001). 

7. Earth-Moon dynamics: The mass of the Earth- 
Moon system, the lunar tidal acceleration, pos- 



sible geocenter variations and related processes 
as well as further effects can be investigated in 
detail. 

8. Time scales: The lunar orbit can also be consid- 
ered as a long-term stable clock so that LLR can 
be used for the independent realization of time 
scales. Those features shall be addressed in the 
future. 




Fig. 7. Distribution of retro-reflectors on the Moon surface. 

To use the full potential of Lunar Laser Ranging, 
the theoretical models as well as the measurements 
require optimization. Using the 3.5 m telescope at 
the new Apollo site in New Mexico, USA, millime- 
ter ranging becomes possible. To allow the determi- 
nation of the various quantities of the LLR solution 
with a total gain of resolution of one order of magni- 
tude, the models have to be up-dated according to the 
IERS conventions 2003, and made compatible with 
the IAU 2000 resolutions. This requires, e.g., to bet- 
ter model 

• higher degrees of the gravity fields of Earth and 
Moon and their couplings; 

• the effect of the asteroids (up to 1000); 

• relativistically consistent torques in the rota- 
tional equations of the Moon; 

• relativistic spin-orbit couplings; 

• torques caused by other planets like Jupiter; 

• the lunar tidal acceleration with more periods 
(diurnal and semi-diurnal); 



• ocean and atmospheric loading by updating the 
corresponding subroutines; 

• nutation using the recommended IAU model; 

• the tidal deformation of Earth and Moon; 

• Moon's interior (e.g. solid inner core) and its 
coupling to the Earth-Moon dynamics. 

Besides modeling, the overall LLR processing shall 
be optimized. The best strategy for the data fitting 
procedure needs to be explored for (highly) corre- 
lated parameters. 

Finally LLR should be prepared for a renaissance 
of lunar missions where transponders (e.g. Degnan 
2002) or new retro-reflectors may be deployed on 
the surface of the Moon which would enable many 
pure SLR stations to observe the Moon. NASA is 
planning to return to the Moon by 2008 with Lu- 
nar Reconnaissance Orbiter (LRO, 2005), and later 
with robotic landers, and then with astronauts in the 
middle of the next decade. The primary focus of 
these planned missions will be lunar exploration and 
preparation for trips to Mars, but they will also pro- 
vide opportunities for science, particularly if new re- 
flectors are placed at more widely separated locations 
than the present configuration (see Fig. 7). New in- 
stallations on the Moon would give stronger determi- 
nations of lunar rotation and tides. New reflectors on 
the Moon would provide additional accurate surface 
positions for cartographic control (Williams et al. 
2005b), would also aid navigation of surface vehicles 
or spacecraft at the Moon, and they also would con- 
tribute significantly to research in fundamental and 
gravitational physics, LLR-derived ephemeris and 
lunar rotation. Moreover in the case of co-location of 
microwave transponders, the connection to the VLBI 
system may become possible which will open a wide 
range of further activities such as frame ties. 

5 Conclusions 

For the IERS, LLR has contributed to the realiza- 
tion of the International Terrestrial Reference Frame 
ITRF2000 and to combined solutions of Earth Orien- 
tation Parameters. 

Additionally, LLR has become a technique for 
measuring a variety of relativistic gravity parameters 
with unsurpassed precision. No definitive violation 
of the predictions from general relativity are found. 
Both the weak and strong forms of the EP are ver- 
ified, while strong empirical limitations on any in- 
verse square law violation, time variation of G, and 
preferred frame effects are also obtained. 



LLR continues as an active program, and it can 
remain as one of the most important tools for testing 
Einstein's general relativity theory of gravitation if 
appropriate observations strategies are adopted and 
if the basic LLR model is further extended and im- 
proved down to the millimeter level of accuracy. 

Additional ranging devices on the Moon would 
have benefits for lunar science, fundamental physics, 
control networks for surface mapping, and naviga- 
tion. Demonstration of active devices would prepare 
the way for very accurate ranging to Mars and other 
solar system bodies. 
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